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Topological entropy is a common measure of the rate of mixing in a flow. It can be 
computed by partition methods, or by estimating the growth rate of material lines or other 
material elements. This requires detailed knowledge of the velocity held, which is not 
always available, such as when we only know a few particle trajectories (ocean float data, 
for example). We propose an alternative approximation to topological entropy, applicable 
to two-dimensional hows, which uses only a hnite number of trajectories as input data. To 
represent these sparse data sets, we use braids , algebraic objects that record how strands, 
i.e., trajectories, exchange positions with respect to a projection axis. Material curves 
advected by the how are represented as simplified loop coordinates. The exponential rate 
at which a braid deforms loops over a hnite time interval as the strands exchange places is 
the Finite-Time Braiding Exponent (FTBE) and serves as a proxy for topological entropy 
of the two-dimensional how. We demonstrate that FTBEs are robust with respect to the 
value of numerical time step, details of braid representation, and choice of initial conditions 
inside the mixing region. We also explore how closely the FTBEs approximate topological 
entropy depending on the number and length of trajectories used. 


In geophysical flows and many other ap¬ 
plications it is important to know where 
things can go, and where they come from. 
This is the study of transport and its 
cousin, mixing. Modern methods are 
most powerful when we know the flow 
perfectly through its velocity field. But 
when the data comes from, say, ocean 
floats, it is very sparse and cannot com¬ 
pletely characterize transport. We use 
tools from the mathematical branch of 
topological dynamics to tease out as much 
information as possible about mixing in 
the flow, in particular the degree of “en¬ 
tanglement” that a set of trajectories 
achieves. We call this measure the Finite- 
Time Braiding Exponent. 
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I. INTRODUCTION 

A common approach to studying kinematics 
of physical fluid flows is through their repre¬ 
sentation as dynamical systems, either as dif¬ 
ferential equations that govern paths of ad¬ 
vected particles, or as maps that connect initial 
and terminal locations of the particles. When 
the dynamics are deterministic, the complexity 
of the flow is connected to the chaotic nature 
of particle trajectories, a phenomenon called 
chaotic advection)r ^ 

In dynamical systems, and especially ergodic 
theory, the amount of complexity present in the 
system is quantified by various notions of en¬ 
tropy. Informally, positive entropy in the sys¬ 
tem is the indicator of sustained complex behav¬ 
ior of trajectories. Formally, there are several 
non-equivalent notions of entropy that can be 
used to characterize a dynamical system. The 
two most common entropies in deterministic dy¬ 
namical systems are metric (Kolmogorov-Sinai) 
entropy and topological entropy. 

Both metric and topological entropies can be 
described by considering a setting in which two 
trajectories of length T can be distinguished 
only if they are further than some resolution e 
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apart. If we count the number of £-distinct tra¬ 
jectories of length T, we expect their number to 
grow as T —> oo. If the growth is asymptoti¬ 
cally exponential ~ e hT , then the rate h is the 
entropy of the flow. The crucial difference be¬ 
tween topological and metric entropy is in the 
way that £-distinct orbits are counted: topolog¬ 
ical entropy counts all e-distinct orbits, while 
metric entropy counts only typical ^-distinct or¬ 
bits, where “typical” allows for ignoring a set of 
orbits that is of zero measure with respect to 
an invariant measure)- 3 This means that each 
invariant measure has an associated metric en¬ 
tropy, and topological entropy is their supre- 
mum. An immediate practical consequence 
is that topological entropy provides an upper 
bound for metric entropy. For a further theo¬ 
retical overview, see ref. |3j for applications, see 
refs. dj and [5] 

Counting distinct orbits to estimate either of 
the entropies is feasible only for certain classes 
of dynamical systemsPHowever, topological en¬ 
tropy can be estimated by discretizing the ad- 
vection operator. This approach has been used 
both for analysis 7 ® and synthesid^ of dynamics. 

Alternatively, metric entropies can be esti¬ 
mated using Lyapunov exponents, which mea¬ 
sure the rate of separation of particles as they 
are advected by the flow. Lyapunov exponents 
are the mean of local deformation rates along 
trajectories of the flow. Local measurements 
of deformation — Lyapunov exponents — and 
global measurements of complexity — entropies 
- are connected by the Margulis-Ruelle in¬ 
equality which states that the spatial mean of 
positive Lyapunov exponents with respect to 
the invariant measure is an upper bound for the 
metric entropy with respect to that measure. 
The Pesin formula strengthens the inequality 
to an equality for those invariant measures that 
are non-singular with respect to the standard 
volume measure.® 10 

In applied problems, trajectories are almost 
always observed over finite times, and thus def¬ 
initions of entropies do not strictly apply. Nev¬ 
ertheless, rates of deformation are used as in¬ 
dicators of complexity in geophysical dynamics, 
in addition to statistical quantities such as rel¬ 


ative dispersion ratef® and effective diffusiv- 
ityP® In particular, the amount of mixing in 
a region is commonly measured by exponen¬ 
tial rates of separation of fluid trajectories, ei¬ 
ther over finite time horizons (Finite-Time Lya¬ 
punov Exponents,^ or FTLEs) or up to finite 
scales (Finite-Size Lyapunov Exponents)-— I— I or 
FSLEs). Computation of FTLE and FSLE 
fields gained prominence in studyingbarriers to 
material transport in fluid flows)®*® as ridges 
and troughs of FTLE/FSLE fields were used as 
proxies®*® to Lagrangian Coherent Structures 
(for an introduction to the topic see ref. IMI) . 
Despite some issues in using FTLE and FSLE 
for detectio n of Lagrangian Coherent Struc- 
tures IAj25|26j these techniques remain popular in 
analysis of geophysical fluid flows P® 

In experimental and numerical fluid flows, 
Lyapunov exponents might be difficult to com¬ 
pute efficiently and reliably as they require 
knowledge of the gradient tensor of the ve¬ 
locity field, and precise temporal resolution of 
data!® As an alternative, topological entropy 
of two-dimensional flows governed by ODEs can 
be estimated by measuring growth of material 
curves, as shown by Newhouse and Pignataro!® 
A material curve is a curve of initial condi¬ 
tions advected as a set by the dynamical sys¬ 
tem. Material curves are expected to grow ex¬ 
ponentially fast in chaotic flows; the largest ex¬ 
ponential rate, over all material curves, is equal 
to the topological entropy in smooth flows, or 
else it provides a lower bound. In practice the 
advection of material curves requires tracking 
an exponentially-growing number of trajecto¬ 
ries, which can exceed a computer’s memory. 

We work at the other extreme of available 
data, assuming that a sparse set of trajectories 
is the only information known about the flow. 
While this regime seems restrictive, such data 
sets are common in physical oceanography, com¬ 
prising trajectories of floats and drifters pas¬ 
sively advected by currents!®- 30 Other exam¬ 
ples include granular m edial crowds of peo¬ 
ple)® and animal flocks)®® all of which can 
behave similarly to particles advected by flows, 
without being driven by physically-observable 
vector fields. Even when detailed models of flow 
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velocities are available, it might be preferable to 
compute coarse estimates based on several tra¬ 
jectories when speed of analysis is important. 
In all these examples, topological approaches 
have proved useful to measure complexity of 
two-dimensional dynamics. 

The topological theory for measuring rates of 
deformation is intuitively similar to direct ma¬ 
terial line advection. However, instead of a de¬ 
tailed bookkeeping of material curves and par¬ 
ticle trajectories, the data is simplified drasti¬ 
cally using two symbolic representations: braids 
and loops. The trajectories are represented as 
braids: time-ordered sequences of symbols that 
encode the manner in which trajectories ex¬ 
change their order along a chosen axis ESEll Ma¬ 
terial lines are represented by topological loops: 
closed “rubber bands” that wrap around trajec¬ 
tories!^ Both braids and loops can be econom¬ 
ically represented using symbols: generators for 
braids, and loop coordinates for loops. Acting 
on loops with braids amounts to studying tight 
material curves that enclose a set of trajectories. 
These curves grow as the trajectories evolve in 
time. Because all this is done using symbolic 
representations, estimates of loop growth can 
be computed very efficiently compared to the 
full model. 

In applied dynamical systems, braids were 
first used to represent dynamical evolution of 
points that lie on periodic orbitsP^®!] When 
trajectories are periodic, recording them over 
one common interval and representing them us¬ 
ing braid generators results in an element of the 
braid group. The entropy of braids can be com¬ 
puted very precisely, sometimes even analyti¬ 
cally (or at least in terms of the largest root 
of a known polynomial). Braids conveniently 
label classes of continuous mappings of two- 
dimensional domains. As topological entropy 
of the flow is associated with its flow map, it 
can be shown that entropy of a braid of a pe¬ 
riodic orbit in such a flow is a lowe r bound for 
topological entropy of the flowP^®H!U 

When measured trajectories are not periodic, 
a true topological representation would involve 
an infinite braid. Nevertheless, we can still esti¬ 
mate exponential rates of stretching of loops us¬ 


ing braids of finite trajectories, similar to calcu¬ 
lations given in ref. [37] We use the name Finite- 
Time Braiding Exponent for the proposed com¬ 
plexity measure, as it intuitively corresponds 
to Finite-Time Lyapunov Exponent calculation. 
However, FTBEs differ from FTLEs in that the 
latter is only a local measure of deformation 
around a single trajectory. 

Computing FTBEs amounts to finding the 
maximal rate of deformation that a loop can 
experience under the action of the braid; how¬ 
ever, not all loops deform at the same rate. 
Allshouse and Thiffeault'Sl showed that slow- 
growing loops often encircle coherent structure 
boundaries. In flows with sparse trajectory 
data, a more detailed understanding of loop de¬ 
formation by braids of trajectories is thus rel¬ 
evant for both identification of coherent struc¬ 
tures and measurement of mixing. 

Other than several graphs in refs. [30 and [37 
there are no results about how FTBEs depend 
on the way in which data is collected, or how 
tightly they might approximate the topological 
entropy of the flow. A theorem in this direc¬ 
tion would be invaluable, but in its absence we 
explore how FTBEs depend on 1. the number 
of trajectories in the data set (Section HID2); 
2. the location of initial conditions for trajecto¬ 
ries; and 3. the length of trajectories (both in 
Section III D 1) . Numerical computations with 
the Aref Blinking Vortex flow indicate that the 
spatial variance of FTBEs decays with both in¬ 
creased number and length of trajectories. The 
mean FTBE increases with the number of tra¬ 
jectories until a saturation is achieved. We 
propose models for this growth in Section |TV} 
as well as compare the saturation bound with 
other estimates of topological entropy. 


II. MATHEMATICAL BACKGROUND 

As mentioned in the introduction, we assume 
that the only data known about a dynamical 
system on a two-dimensional domain is a finite 
number of trajectories of equal, finite lengths. 
In this section, we describe how to represent 
the trajectories as braids and how to calcu- 
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late Finite-Time Braiding Exponents from the 
braid. Finally, we briefly describe braidlab, 4 ' 
a publicly-available MATLAB toolbox imple¬ 
menting the calculation of braids and FTBEs. 


A. Representation of braids and loops 

Our dynamics will be assumed to be defined 
on a closed unit disk D in a plane. Since the 
braid theory approach uses only topological ar¬ 
guments, any other domain homeomorphic to a 
disk, i.e., that can be continuously deformed to 
a disk, will require no special treatment. 

A physical braid (or geometric braid) is a col¬ 
lection of finitely-many planar continuous tra¬ 
jectories on a (bounded) time interval / = 
[to, to + T]. When they are thought of as curves 
in the extended three-dimensional space D x /, 
trajectories are also called strands (or strings). 
At any time t € J, strands evaluate to a set of 
n distinct points, often called punctures , which 
move around in time. 

Constructing a braid from a physical braid is 
a method of discarding information about ge¬ 
ometry (e.g., distances between strands), while 
retaining information about topology (e.g., rel¬ 
ative locations of the strands). We keep only 
those time instances, called crossings , in which 
there is a substantial change in the positions of 
strands. Crossings are defined with respect to 
a projection axis — a line passing through the 
origin of the domain at an angle a\ 36 l 37 l Strands 
are ordered and indexed 1,2 ,,n according to 
their projection onto the axis. As the punctures 
move in time, they will exchange their indices 
when they cross. 

A braid b is a sequence of symbols b = 
t>i &2 • • • where each symbol corresponds to a 
crossing in a physical braid. The symbols are 
taken from the set of 2 (n — 1) braid generators 
{erf 1 , af 1 ,..., cr^ i 1 _ 1 ^}. Each symbol erf 1 cor¬ 
responds to a crossing of strands i and i + 1, 
either clockwise (+1) or anticlockwise (—1).® 

Loops are counterparts of advected material 
lines, in the same way as braids are counter¬ 
parts of trajectories. To specify a loop in the 
disk, draw one or more closed curves that do 


not intersect themselves, each other, the bound¬ 
ary, nor pass through one of the punctures. In 
the topological picture, curves can be continu¬ 
ously deformed into each other without passing 
through one of the punctures. This allows us 
to “tighten” each loop around the punctures it 
encloses. Loops that can be tightened to a sin¬ 
gle puncture or expanded to the disk boundary 
are never affected by the motion of punctures; 
we do not consider those loops. The remaining 
( essential ) loops are deformed by the motion as 
if they were rubber bands caught on the punc¬ 
tures. 

We encode loops with a particular integer co¬ 
ordinate system, introduced by Dynnikov 48 and 
explained in detail in refs.[371 and (49l Dynnikov 
coordinates provide a one-to-one and onto cor¬ 
respondence between essential loops and vec¬ 
tors in Z 2n_4 . Each braid generator af 1 acts 
on loops via a piecewise-linear function erf 1 : 
Z 2n — 4 —>. Z 2 " -4 , while the action of full braid is 
formed by composition of such functions. The 
functions are given in ref. [37] and 33 

Using braids and Dynnikov coordinates we 
represent the dynamics of material lines by 
integer-valued functions. Given a braid b and 
a vector v € Z 2ra-4 encoding a loop, we will 
use multiplication b • v (or br>) to denote the 
nonlinear action of b on v. 


B. Finite-Time Braiding Exponents 

The word length L(h) of a braid b = 
b\b 2 ■ ■ ■ bk is k. While trivial to compute,® the 
word length does not capture topological com¬ 
plexity very well. For example, imagine a set of 
strands that are sparsely placed on the domain, 
and dynamics that only “wiggles” them slightly. 
With certain choices of projection angle a, wig¬ 
gles will produce plenty of crossings among the 
strands^ with no topological complexity what¬ 
soever.® 

Braid entropy is a measure of braid complex¬ 
ity derived from topological entropy of continu¬ 
ous maps, specifically pseudo-Anosov maps!- 51 1 
In topological surface dynamics, braids are 
used as labels for classes of continuous maps 
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(a) e E 



(b) (TlfTl ■ t E 


FIG. 1. Top: The loop £ E for the 3-punctured disk, 
used in the iterative computation of topological en¬ 
tropy 0 - The three main punctures are black, 
while the white puncture serves as the “basepoint” 
for the loop. Bottom: i E after deformation by the 
braid cricri, which exchanges the first pair of punc¬ 
tures twice, clockwise. 


that can be continuously deformed into each 
other (homotopy-equivalent homeomorphisms, 
see refs. [38] [52] and [53]) . The topological en¬ 
tropy of a homeomorphism is directly connected 
to the maximal asymptotic exponential rate of 
stretching that a curve on a surface will expe¬ 
rience un der r epeated action of the homeomor- 
phismP ^ 9 ! 51 * 54 !^ The topological entropy of a 
braid h(h) is then the minimal entropy of a 
homeomorphism out of the entire class of home¬ 
omorphisms labeled by the braid b. 

ThurstoiPS showed that h{ b) can be com¬ 
puted from the growth rate of a loop £e under 
repeated action of the braid b: 


M b ) = 


1. |b He 

, hm 7 lo SVr 

k —^oo /c \^E\ 


(i) 


Here we take the loop 4 to be a generating set 
for the fundamental group of the n-punctured 
disk OP^It can be represented in Dynnikov co¬ 
ordinates by adding an extra “basepoint” punc¬ 
ture that does not participate in the braid b 
(Figure 1). The length |-1 is the number of in¬ 
tersections of the loop with the horizontal axis 


to the left of the basepoint (dashed line in Fig¬ 


ure 1). 


This iterative calculation may however be 
an inappropriate measure of complexity of tra¬ 
jectories. First, unless the trajectories corre¬ 
spond to periodic orbits, there is no justifi¬ 
cation for repeated application of the braid, 


since in that case the physical braid does not 
represent a homeomorphism of a fixed punc¬ 
tured surfaced 7 Second, flows that stretch ma¬ 
terial sub-exponentially will always be assigned 
h(h) = 0, which may be too crude of a measure 
in practice, e.g., for estimating time-scales for 
transport. 

Dynnikov and WiespSI used a non-iterated 
version of 0 to study the connection between 
algebraic and geometric properties of braids. 
We adapt their definition here to define the 
Finite-Time Braiding Exponent as a measure 
of complexity of trajectories. 

Definition (Finite-Time Braiding Exponent). 
Let b be a braid corresponding to n trajectories 
over a time interval of length T. The Finite- 
Time Braiding Exponent is given by 

FTBE(b) = 1 log ^ (2) 

where £e is the loop representing a generating 
set of the non-oriented fundamental group on 
n-punctured disk. 

A similar definition of exponential deforma¬ 
tion rate appeared in ref. 1371 which uses an ar¬ 
bitrary loop £ instead of £e- 


III. CONSTRUCTING BRAIDS AND 
COMPUTING FTBES 

Our ultimate goal is to use FTBEs as an ap¬ 
proximation to topological entropy of a planar 
dynamical system. In this section we discuss 
the issues involved in constructing braids from 
data, introduce our model system, and examine 
the sensitivity of braids and FTBEs to various 
model parameters. 


A. Practical considerations 

To compute FTBEs, we must first construct 
braids. To do this, we record a finite number of 
concurrent, non-intersecting trajectories of the 
dynamics. Such a data set is determined by 
three parameters: 
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1) the number n of trajectories, 

2) the initial conditions of trajectories, and 

3) the time interval [to, to+T\ during which tra¬ 
jectories are recorded. 

The MATLAB toolbox braidlatP^ converts 
trajectories into braids and then computes 
FTBEs. braidlab is an open-source library 
that implements a wide set of operations on 
braids, loops, and related structures. In par¬ 
ticular, braids generated from data as discussed 
here are stored as databraid data structures, 
which record both crossings of strands as gener¬ 
ators and the times at which crossings occurred. 

The computational time to convert n trajec¬ 
tories to a braid scales quadratically with n, and 
linearly with both the length of trajectories and 
the total number of crossings between strands. 
The algorithm depends on two additional pa¬ 
rameters: 

4) the angle a of the projection axis used to 
assign order to strands; and 

5) the time step r at which trajectories are sam¬ 
pled, assuming uniform sampling. 


t t + T 

(a) Identity (trivial) braid 



(c) <Ji(Ti braid 

FIG. 2. Three physical braids with same orders of 
strands i, i + 1 at sampling times t and t + r, but 
different sequences of generators. Using time step 
t/2 would resolve the crossings correctly. 


B. The Aref Blinking Vortex flow 


To accurately convert discretized trajectories to 
a braid, braidlab requires that strands cross 
only between sampled times, and that any two 
strands cross at most once per time step. If 
strands cross exactly at the moment of sam¬ 
pling, a braid generator cannot be unambigu¬ 
ously assigned. Such degeneracies may be re¬ 
solved by varying the angle a.® 

When two strands cross more than once in 
a single time step, other ambiguities can arise. 
For example, |Figure 2| shows three different 
physical braids with same strand positions at 
the sampling time steps. The ambiguity can 
clearly be resolved by reducing r. 

The sections that follow demonstrate the de¬ 
pendence of FTBEs on each of the five param¬ 
eters listed above: number of strands n, initial 
conditions, duration of time interval T, projec¬ 
tion angle a, and time step r. The input tra¬ 
jectories are generated by a mixing dynamical 
system, the Aref Blinking Vortex flow. 


The Aref Blinking Vortex flow® is an ideal¬ 
ization of a device that stirs fluid using two 
counter-rotating vortices in a circular domain. 
The vortices are otherwise identical, and are 
positioned along the diameter of the unit cir¬ 
cle at distances ±b from its center. They are 
turned on alternately (“blinked”) during half 
of the period Tp of the protocol, resulting in 
continuous, piecewise-differentiable trajectories 
of fluid parcels, shown in Figure 3 After fix¬ 
ing the geometry, and period of the protocol 
to Tp = 1, the only free parameter is the non- 
dimensional circulation fi. We restrict our anal¬ 
ysis to /r € [3,9] for which numerical experi¬ 
ments indicate that the flow is mixing in the 
entire domain. 1 


We quantify the complexity of the flow by 
estimating topological entropy h using direct 
advection of material curves“ Numerically, we 
seed initial conditions along a randomly selected 
straight material line and advect those points 
forward. When distances between neighboring 
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FIG. 3. Material advection patterns and two parti¬ 
cle trajectories of the Aref Blinking Vortex flow in 
the mixing regime. 


points grow too large, we linearly interpolate 
additional points between them to ensure that 
details of bends in the material curve are well 
represented. 

|Figure 4| shows the estimated values of topo¬ 
logical entropy h that will be used as reference 
for our later calculations of FTBEs. Topological 
entropy h is approximately linear in circulation 
/i, which determines the dominant time scale of 
the flow. 



Circulation /x 

FIG. 4. Topological entropy of counter-rotating 
Aref Blinking Vortex flow 


C. Robustness of braid construction 


1. Time step 


Let us investigate how the choice of the sam¬ 
pling interval r affects accuracy of constructed 
braids and their FTBEs. We will use the word 
length of braids as a proxy for the accuracy of 
constructed braids, as reducing r increases the 
number of crossings in the braid. 

Vortices in Aref Blinking Vortex flow are 
modeled by singularities which we slightly reg¬ 
ularize to avoid infinite velocities. After regu¬ 
larization, time step r* = 10~ 4 is sufficient to 
accurately sample trajectories of the flow in the 
studied regime of circulations. For each set of 
initial conditions we compute a braid 6* with 
time step r*, which serves as a reference, and 
additional braids b T with larger values of time 
step r. Figure [5] shows how the relative errors 



FTBE(b r ) 

FTBE(b*) 


, and 

, L( b r ) 


L( b*) 


(3) 


depend on time step r. To reduce the effect of 
spatial variations, we plot the largest relative er¬ 
rors with among different uniformly-initialized 
braids. 



FIG. 5. Effect of time step t on relative error in 
FTBE and braid word length pi with respect to 
the reference time step t* = 10“ . Graph shows the 
maximum error among 100 braids of 50 uniformly 
initialized strands. 


In this section we investigate numerical issues 
involved in braid construction, which could lead 
to poorly-defined braids if not treated properly. 
The two main issues are the choice of time step 
(sampling rate, section III C 1) and the choice 
of projection angle (Isection III C 21) . 


Figure [5] suggests that FTBEs are more re¬ 
silient to large time steps than word length. 
This implies that some parts of the braid do 
not contribute significantly to deformation of 
loops. For example, a braid of four genera¬ 
tors CTiCr^ViiT” 1 introduces no net deformation 
as each inverse generator reverses the action of 




















the previous one. Such sequences can be gener¬ 
ated by pairs of strands that are neighbors along 
the projection line, but are otherwise located far 
from each other in the domain. 

As expected, the effect of large time steps is 
more pronounced for flows with larger circula¬ 
tions /it, as particles move farther within each 
time step. We use r = ICE 2 for the rest of the 
paper, which is appropriate for our range of /i. 


2. Angle of the projection axis 


The angle a determines the projection line 
used to assign order indices to strands of a 
physical braid. As mentioned at the opening 
of Section m a may need to be adjusted to 
ensure that trajectories can be converted to a 
braid. We show that even though the braid 
itself changes depending on a, the effect on 
FTBEs is minimal. 


When a physical braid consists of periodic 
trajectories, changing a results in conjugation 
of the corresponding braid b by another braid 
c a . In other words, dependence on angle a en¬ 
ters into the braid as c” 1 bc a . However, when 
strands in the physical braid are not periodic, 
changing a is not guaranteed to affect the braid 
through conjugation onlyP-^ Nevertheless, since 
changing a only rotates the coordinate sys¬ 
tem, the FTBEs should not significantly depend 
on a. 


Figure 6 shows how length T of trajectories 


influences statistics of FTBEs and braid word 
length with respect to projection angle a. It is 
clear that the change of projection angle min¬ 
imally affects the value of FTBEs for longer 
braids. Over the entire range of times T simu¬ 
lated, the Relative Standard DeviatioiP^ (RSD) 
of FTBEs decays inverse-proportionally with 
length of the physical braids, while the standard 
deviation of length remains approximately con¬ 
stant, indicating that FTBEs are more robust 
than braid length with respect to change of a. 
Based on this, we set a = 0 in the remainder of 
this paper. 



(a) Mean FTBE 



10 ° 10 1 10 2 
Trajectory length T 

(b) Relative Standard Deviation of 
FTBE 



(c) Mean braid word length 


O n = 3 

0 n = 25 

v n = m ( 

-0 slope 
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00 rfp 
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10 ° 10 1 10 2 
Trajectory length T 


(d) Relative Standard Deviation of 
braid word length 


FIG. 6. Effect of projection angle a on FTBEs and 
braid word length of a single physical braid. Statis¬ 
tics are computed with respect to uniform distribu¬ 
tion of 100 angles a G [0,7r). Circulation is [i = 5. 
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D. Parameters for braid construction 


In lsection III Cl we addressed numerical issues 
that can lead to deviations from the idealized 
mathematical representation of braids. In this 
section we turn to parameters that affect the 
nature of the braid we obtain. These are the 
choice of initial conditions, the length (dura¬ 
tion) of trajectories (section III Dll, and the 
number of trajectories (section IIID 2). 


1. Initial conditions and length of trajectories 


FTBEs can be interpreted as measures of av¬ 
erage deformation of a loop over finite segments 
of physical braids. Averaged measurements in 
chaotic systems quickly stabilize to constant 
values independent of precise initial conditions, 
which is formally asserted by various flavors of 
ergodic theorems. Since numerical evidence in¬ 
dicates that the Aref Blinking Vortex flow is 
mixing for /i £ [3,9], the FTBEs should not de¬ 
pend on initial conditions for long-enough tra¬ 
jectories. 

|Figure 7| confirms the expected ergodic be¬ 
havior. After the rapid initial transient, the 
spatial mean of FTBEs settles to a constant 
value, while spatial fluctuations decay following 
a power law. Since the length of the transient 
does not depend on number n of strands, we 
set the length of trajectory to T = 100 in fu¬ 
ture simulations. 

The slope —1/2 of decay of fluctuations is 
commonly associated with the Central Limit 
Theorem for sums of random variables. Vari¬ 
ants of the Central Limit Theorem apply to the 
decay in fluctuations of Fi nite-T ime Lyapunov 
Exponents in mixing flows.® 61 While we pro¬ 
vide no rigorous proof that the same applies 
to FTBEs, it is not surprising to see the same 
behavior since FTBEs and FTLEs are closely- 
related quantities. 



(a) Mean FTBE 



(b) Relative Standard Deviation of 
FTBE 

FIG. 7. Effect of trajectory length T on FTBE for 
an ensemble of physical braids. Statistics are com¬ 
puted with respect to uniform initialization of 100 
physical braids. The circulation is g = 5. 


2. Number of strands 


It has been observed that using more 
traje ctories increases the complexity of the 
braid.®®® To study how FTBEs depend on 
the number of strands n, we fix the length of 
trajectories to T — 100 and at each value n 
compute mean and standard deviation 59 - of the 
FTBE with respect to initial conditions. 


Figure 8 shows that mean FTBE and mean 


braid word length grow in a different manner 
with addition of strands. The mean FTBE ap¬ 
pears to limit to topological entropy h; on the 
other hand, the number of crossings L is un¬ 
bounded, growing proportionally to the number 
of pairs of strands n 2 . The standard deviations 
of both FTBE and word length decay according 
to a power law, although with slightly different 
slopes. 

To generate |Figure 8 we computed approxi¬ 
mately S x ^points n « 150, 000 trajectories of 
the Aref Blinking Vortex flow. While such a 





























10 



(a) Mean FTBE. Horizontal lines 
indicate topological entropies h. 



(b) Relative Standard Deviation of 
FTBE 




(d) Relative Standard Deviation of 
braid length 


available experimentally is often much smaller, 
especially in systems featuring human crowds, 
flocks of birds, or oceanic floats. We can, how¬ 
ever, generate collections of braids by re-using 
data from a much smaller set of trajectories. 

If we have access to m measured trajectories, 
we can form a single m-stranded braid, but also 
select subsets of n < to trajectories to form 
smaller n-subbraids. There are different 
n-stranded subbraids, which is a number that 
grows fast as n —¥ m/2 from above or below. 
Therefore, it is easy to obtain large collections 
of n-braids by selectively including strands of a 
bigger set of to trajectories. 

Despite the practical appeal, it is not clear 
that statistics of FTBEs and braid lengths will 
be similar if calculated using either indepen¬ 
dent braids or subbraids of a larger braid; sub¬ 
braids are not mutually independent, as they 
re-use the same set of trajectories from the orig¬ 
inal data set. To compare statistics, we gen¬ 
erate subbraids from a single set of to = 550 
trajectories for the same values of n used in 
the calculations |Figure 8 which used indepen¬ 
dent sampling. Figure |9] overlays the results 
of two calculations and indicates that statis¬ 
tics of FTBEs over subbraids match almost per¬ 
fectly with statistics over independently sam¬ 
pled braids. 

An additional practical benefit of using sub¬ 
braids is the speed of generating them. The 
time to extract one n-subbraid from a full m- 
braid scales with the braid word length of the 
full braid. This procedure is generally much 
faster than forming an independent n-braid, 
even if the time to simulate new trajectories is 
negligible. In our example statistics over sub¬ 
braids were computed 5 to 10 times faster than 
over independent braids. 


FIG. 8. Influence of number n of strands on FTBE 
and length of an ensemble of physical braids. Statis¬ 
tics are computed with respect to uniform initial¬ 
ization of 100 physical braids. 


large number of trajectories can be computed 
for numerical flows, the number of trajectories 


IV. EXTRAPOLATING TOPOLOGICAL 
ENTROPY FROM FTBES 

|Figure 10| demonstrates that topological en¬ 
tropy h of the flow and mean values of FTBEs 
are strongly correlated, therefore the FTBE can 
serve as a finite-scale proxy for h. To go further, 
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(a) Mean FTBE 



(b) Relative Standard Deviation of 
FTBE 



FTBE 


(c) Distribution of FTBE at 
n = 500, value with largest 
difference in RSD between two sets. 


no formal verifications.^ The results of Sec¬ 
tion EE D 2| suggest that FTBEs behave simi¬ 
larly. 



FIG. 10. Correlation of topological entropy and 
mean FTBEs computed using different numbers of 
strands n. Points were obtained by varying circula¬ 
tion /r. 


To investigate the approach of FTBEs to 
topological entropy h, we compute the relative 
difference 1 — FTBE / h, and compute its statis¬ 
tics across sets of S = 100 subbraids of a phys¬ 
ical braid of to = 550 trajectories. |Figure lT] 
shows that the relative difference approximately 
decays as a power law with slope —1/2 for a 
large range of strands n. Nevertheless, the de¬ 
cay deviates from the power law as n increases, 
which means that extrapolating data using a 
power law model would overestimate the value 
of topological entropy. 


FIG. 9. Comparison of FTBE for braids sampled 
independently and braids formed as subbraids of m 
strands chosen from a single m = 550 braid. Sample 
size at every parameter value is S = 100, circulation 
H = 9. 


and extrapolate values of FTBE to obtain the 
actual value of topological entropy h, we need 
that:l) the FTBE indeed approaches topologi¬ 
cal entropy as n —>• oo, and 2) the growth of 
FTBE can be described by a relatively simple 
functional model, whose parameters can be es¬ 
timated. Neither of these statements are rigor¬ 
ously known to be true. Previous numerical re¬ 
sults indicate that topological entropy of braids 
formed by periodic trajectories approaches the 
topological entropy of the flow, but there are 



FIG. 11. Relative difference of the mean FTBE and 
topological entropy h for the Aref Blinking Vortex 
flow. Statistics are computed with respect to 100 
n-subbraids of a set of 550 trajectories of length 
T = 100. 

To estimate the decay of the relative differ¬ 
ence more carefully, we fit two functions to the 
calculated points: 
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Tempered Power Law: 


n?e Pn 


(4) 


Log-normal CDF: 

erfc{ 7 log fln} 


(5) 


These two functions are often used as compet¬ 
ing statistical models for tails of distributions, 
where the criterion for choosing one over the 
oth er is b ased on a maximum likelihood estima- 
tionPSUl Since relative difference at n strands 
is not associated with a distribution, we use the 
more classical weighted least-squared-error fit¬ 
ting of models to data. Weights assigned to 
data points are inversely proportional to the 
standard deviation of FTBEs, which enforces a 
tighter fit to points at larger n, due t o dec aying 
standard deviation, shown in Figure [9(b) [ 
Figure [12] shows the quality of fit of the two 
models to values of FTBEs. Both functions fit 


well to dat a, as evidenced by Figure 12(a) Fig¬ 
ure |12(b)| shows residual errors for both mod¬ 
els at different values of circulation and demon¬ 
strates that the tempered power law fits slightly 
better. 

The true test of the models should be their 
ability to correctly predict the topological en¬ 
tropy h as their limit. Instead of using values of 
h as a known limit of proposed models, we intro¬ 
duce a limit parameter and allow it to vary 
during fitting. This additionally means that we 
fit our model functions to FTBE data instead 
of the relative difference to a reference value. 
Other than these modifications, models and fit¬ 
ting technique remain the same. 

Unsurprisingly, we again find that both mod¬ 
els provide a tight fit to data; however, their 
limits differ due to different rates of growth as 
n —> oc. |Figure 13] shows relative differences 
1 — hoo/h between extrapolated limits of FTBE 
hoc and topological entropy h, depending on cir¬ 
culation /r. In addition to fitting models to all 
available data points, we also fit them to data 
sets with points at higher values of n removed, 
to test robustness of the extrapolations hoo to 
the amount of data used. Only non-negative 



io‘ io 2 

n strands 


(a) Fit of functions ^ and 
(lines) and data (circles) to relative 
difference between mean FTBE and 
topological entropy 
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(b) Weighted least-squared-error 
fit residuals 


FIG. 12. Quality of fit of Tempered Power Law and 
Log-normal CDF models to computed mean FTBE 


values of differences are consistent with the ex¬ 
pectation that values of FTBE approach h from 
below as n —> oo. 

Figure 13| demonstrates that the Log-normal 
model consistently overestimates the value of 
topological entropy. On the other hand, extrap¬ 
olations of the tempered power law provide a 
better estimate than the raw FTBE calculation 
at the largest number of strands (in blue), while 
remaining consistent with the theory. We note 
that extrapolating the fit up tom = 350 strands 
yields an estimate that is of similar quality to es¬ 
timation of h, by FTBE at n = 500 strands with¬ 
out extrapolation. This translates to savings 
in data collection and analysis time, since con¬ 
verting trajectories to braids takes 0(n 2 ) time 
to complete, and computing a single braid at a 
large n can take longer than computing many 
collections of subbraids at m < n. 

The numerical evidence indicates that a tem¬ 
pered power law models the approach of FTBE 
to topological entropy, with power law expo- 
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(a) Tempered Power Law 



FIG. 13. Comparison of h obtained by extrapo¬ 
lating tempered power law and log-normal models 
for mean FTBE. Blue lines show the best estimate 
of topological entropy obtained without extrapola¬ 
tion. Different gray lines connect values hoo ob¬ 
tained from models fit to larger and larger point 
sets; number next to the line is the number of 
strands of last point used for fitting. 


nents 7 estimated to lie in [—0.4,—0.5]. Fix¬ 
ing the exponent of the tempered power law 
to 7 = — 1/2 minimally changes quality of fit 
and extrapolation by the tempered power law 
model, suggesting that growth of FTBE follows 
an universal behavior, regardless of circulation 
of the flow. 

We caution, however, against making 
stronger arguments based on the experiments 
performed, because the reference value for 
topological entropy h is an estimate itself, 
computed by material advection. Since the 
precision achieved by FTBE is within 1% to 
5% of h, it is likely that it approaches the con¬ 
fidence interval for the estimation by material 
advection. For quantitative conclusions about 
the approach of FTBE to topological entropy 
of the flow, analytical justification for a model 
of growth of FTBE is needed. 


V. DISCUSSION 

We briefly summarize our findings about 
Finite-Time Braiding Exponents. 

• FTBEs are very robust with respect to the 
projection angle a used to construct the 
braids. 

• Reasonably-small values of the time step 
r lead to well-defined braids, but FTBEs 
may tolerate a larger r. 

• Increasing the length of trajectories in 
mixing dynamics results in convergence 
of FTBEs. At the same time, the spa¬ 
tial mean of FTBEs stabilizes, while their 
spatial fluctuations decay. 

• The spatial mean of FTBEs correlates 
strongly with the topological entropy of 
the flow; as the number of strands n is 
increased, the mean approaches the topo¬ 
logical entropy. 

The Finite-Time Braiding Exponent of n tra¬ 
jectories sampled from the flow can be inter¬ 
preted as a “multi-scale” measure of complexity 
of the flow at scales determined by n. An im¬ 
portant finding is that to estimate complexity at 
different scales, a single collection of n trajecto¬ 
ries can be re-used by sampling subbraids. Even 
though such subbraids are not fully indepen¬ 
dent, since a given trajectory will participate 
in many braids, the procedure yields statistics 
of FTBEs similar in quality to those computed 
using independently-sampled braids. This tech¬ 
nique greatly reduces the required amount of in¬ 
put data without sacrificing the quality of com¬ 
puted FTBEs at various scales. 

The value of FTBEs seems to approach the 
topological entropy of the flow as n is increased, 
but more work is needed to establish the asymp¬ 
totics of that approach. Fitting plausible func¬ 
tional models suggests that for trajectories of 
mixing dynamical systems the difference be¬ 
tween FTBEs and topological entropy decays 
according to a power law tempered by an expo¬ 
nential. If confirmed analytically, this relation 
would allow one to estimate free parameters of 
the model from values of FTBEs at several dif¬ 
ferent scales n, and extrapolate the model as 
n —► 00 to approximate the topological entropy 
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very closely. 

While these conclusions are based on the 
study of the Aref Blinking Vortex flow, it is 
reasonable to expect similar results for other 
mixing flows. In flows that are not spatially 
homogeneous, e.g., those containing coherent 
substructures, lengths of strands and locations 
of their initial conditions will be of crucial im¬ 
portance for interpretation of FTBE values. In 
such situations, other quantities based on braids 
can be computed to provide a more complete 
picture of dynamics. In particular, a promising 
approach may be to couple braid approxima¬ 
tions of coherent structured^ with FTBE esti¬ 
mates of complexity inside those structures to 
assess the spatial distribution of different dy¬ 
namics in more complicated flows, even when 
the only information available is a finite set of 
trajectories. 
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